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Abstract 

I review recent worlc and some new results, performed in collaboration with G. 
Sierra, on the Real-Space Renormalization group method applied to quantum spin 
lattice systems mainly in spatial dimensions one and two, and to spin ladders which 
are somehow in between. The first part of these notes is devoted to non-interacting 
systems in ID and 2D and the role played by the correlations between blocks. The 
second part comprises interacting systems in ID, spin ladders and 2D using the 
standard BRG method. 



Proceedings of the El Esconal Summer School 1996 on STRONGLY CORRELATED 
MAGNETIC AND SUPERCONDUCTING SYSTEMS 

I. INTRODUCTION: BRIEF HISTORY OF REAL-SPACE RG METHODS 

The Real-Space Renormalization Group Method has undergone a great revival since 
the Densty Matrix RG was introduce by White in 1992 [|l|, ||2|. Nowdays, this method is 
considered a powerful numerical tool to get non-perturbative results, specially for interacting 
systems in ID although some recent advances for 2D systems have been obtained [Q, [0- 
Despite being numerical, the DMRG has also been a source of inspiration for analytical 
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studies and this is the framework of the present notes. Moreover, much of the El Escorial 
Summer School has been devoted to real-space RG methods 0, |Q, 0. 

The Renormalization Group method has become one of the basic concepts in Physics, 
ranging from areas such as Quantum Field Theory and Statistical Mechanics to Condensed 
Matter Physics. The many interesting and relevant models encountered in these fields are 
usually not exactly solvable except for some privileged cases in one dimension. It is then 
when we resort to the RG method to retrieve the essential features of those systems in 
order to have a qualitative understanding of what the physics of the model is all about. 
This understanding is usually recasted in the form of a RG-fiow diagram were the different 
possible behaviours of the model leap to the eyes. 

Many authors in the past have contributed significantly to the idea of renormalization 
and it is out of the scope of these notes to give a detailed account on this issue. 

We shall be dealing with the the version of the RG as introduced by Wilson and 
Anderson in their treatment of the Kondo problem, and subsequent developments of 
these ideas carried out by Drell et al. at the SLAG group and Pfeuty et al. [|TU|. 

Real space Renormalization Group (RG) methods originated from the study of the Kondo 
problem by Wilson [Q. It was clear from the beginning that one could not hope to achieve 
the accuracy Wilson obtained for the Kondo problem when dealing with more complicated 
many-body quantum Hamiltonians. The key difference is that in the Kondo model there 
exists a recursion relation for Hamiltonians at each step of the RG-elimination of degrees of 
freedom. The existence of such recursion relation facilitates enormously the work, but as it 
happens it is specific of impurity problems. 

From the numerical point of view, the Block Renormalization Group (BRG) procedure 
proved to be not fully reliable in the past particularly in comparison with other numerical 
approaches, such as the Quantum MonteGarlo method which were being developed at the 
same time. This was one of the reasons why the BRG methods remained undeveloped 
during the '80's until the begining of the '90's when they are making a comeback as one 
of the most powerful numerical tools when dealing with zero temperature properties of 
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many-body systems, a situation where the Quantum MonteCarlo methods happen to be 
particularly badly behaved as far as fermionic systems is concerned [Q. 

As it happens, the BRG gives a good qualitative picture of many properties exhibited 
by quantum lattice Hamiltonians: Fixed points, RG-flow, phases of the system etc. as well 
as good quantitative results for some properties such as ground state energy and others 0, 
li, [ini], |20|. However in some important instances the BRG method is off the correct 



values of critical exponents by a sensible amount. 

The origin of the density matrix RG method relies on the special treatment carried out 
by White and Noack [|l^ on the ID tight-binding model, the lattice version of a single par- 



ticle in a box. It was Wilson [jT8[ the first to point out the relevance of this simple model 
in understanding the sometimes bad numerical performance of the standard Block Renor- 
malization Group (BRG) method. In reference the authors proposed a method called 
Combination of Boundary Conditions (CBC) which performs extremely well as compared 
to the exact known solution of the model. Recently, we have studied the role played by the 
boundary conditions in the real-space renormalization group method by constructing a 
new analytical BRG-method which is able to give the exact ground state of the model and 
the correct 1/A^^-law for the energy of the first excited state in the large iV(size)-limit. 
Yet another branch of applications of DMRG inspired ideas is to use the Superblock for- 
malism Ijl^ without resorting to a Density Matrix. For instance, in it has been found 



that the application of this formalism to the anisotropic Heisenberg model in ID successfully 
improves the standard BRG results of Rabin ||2^ . 



The Density Matrix RG has been originally devised to deal with quantum lattice Hamil- 
tonians. However, recent new applications have been developed by Nishino and coworkers 



|2^, to address the renormalization of classical lattice models. One of the outcomes 
of these studies has been to state the relationship between Baxter's corner transfer matrix 
formalism and the Density Matrix RG of White (in ID). 

The number of new developments on this subject is constantly growing and it is not 
possible to give a full account of all of them here. More applications can be found in the 



rest of contributions to these proceedings devoted to real-space RG methods. 



II. REVIEW OF STANDARD BLOCK RENORMALIZATION GROUP 

METHODS (BRG) 

For the sake of completeness and to set up the notation used throughout these notes, the 
block renormalization group method is revisited in this section along the lines of a new and 
unified reformulation of it based on the idea of the intertwiner operator T to be discussed 
below. This treatment is by all means equivalent to the standard approach presented by 
S.R. White in his contribution to this volume. This formulation has recently allowed us to 
introduce the new Variational and Fokker-Planck DMRG methods il^ on equal footing as 



the standard BRG method. For a more extensive account on this method we refer to [T^ 



and chapter 11 of reference and references therein. 

Let us first summarize the main features of the real-space RG. The problem that one 
faces generically is that of diagonalizing a quantum lattice Hamiltonian H, i.e., 

H\^>=E\^> (1) 

where {ip > is a state in the Hilbert space Ti. If the lattice has N sites and there are k 
possible states per site then the dimension of H is simply 

dimH = (2) 

As a matter of illustration we cite the following examples: /c = 4 (Hubbard model), /c = 3 
(t-J model), k = 2 (Heisenberg model) etc. 

When N is large enough the eigenvalue problem (0) is out of the capability of any human 
or computer means unless the model turns out to be integrable which only happens in some 
instances in ci = 1. 

These facts open the door to a variety of approximate methods among which the RG- 
approach, specially when combined with other techniques (e.g. numerical, variational etc.). 
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is one of the most relevant. The main idea of the RG-method is the mode ehmination or 
thinning of the degrees of freedom followed by an iteration which reduces the number of 
variables step by step until a more manageable situation is reached. These intuitive ideas 
give rise to a well defined mathematical description of the RG-approach to the low lying 
spectrum of quantum lattice hamiltonians. 

To carry out the RG-program it will be useful to introduce the following objects: 

• H, : Hilbert space of the original problem. 

• Ti!: Hilbert space of the effective degrees of freedom. 

• H: Hamiltonian acting in Ti. 

• H': Hamiltonian acting in Ti! (effective Hamiltonian). 

• T : embedding operator : H' — H 

• :truncation operator : H — H' 

The problem now is to relate H, H' and T. The criterium to accomphsh this task is 
that H and H' have in common their low lying spectrum. An exact implementation of this 
is given by the following equation: 

HT = TH' (3) 

which imply that if "i/'^, is an eigenstate of H' then T^'^;, is an eigenstate of H with the 
same eigenvalue (unless it belongs to the kernel of T: T^'^, — 0), indeed, 

HT-^'^, = TH'^i'e, = E'T-^'e, (4) 

To avoid the possibility that T^' — with ^' ^ 0, we shall impose on T the condition, 

T^T = 1h' (5) 

such that 
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^ = T^' ^ ^' = (6) 

Condition (|5|) thus stablishes a one to one relation between Ti' and Im(T) in Ti. 
Observe that Eq. (||) is nothing but the commutativity of the following diagram: 

w ^ n 

H' i i H 

w ^ n 

Eqs. and characterize what may be called exact renormalization group method 
(ERG) in the sense that the whole spectrum of H' is mapped onto a part (usually the bottom 
part) of the spectrum of H. In practical cases though the exact solution of Eqs. (^ and (H) 
is not possible so that one has to resort to approximations (see later on). Considering Eqs. 
(0) and (^) we can set up the effective Hamiltonian H' as: 

H' = T^HT (7) 

This equation does not imply that the eigenvectors of H' are mapped onto eigenvectors of 
H. Notice that Eq.(^ together with (|^) does not imply Eq. (^). This happens because 
the converse of Eq.(|^), namely TT"!" ^ l-^ is not true, since otherwise this equation together 
with (^) would imply that the Hilbert spaces Ti and Ti' are isomorphic while on the other 
hand the truncation inherent to the RG method assumes that dimTi' < dimTi. 

What Eq. (0) really implies is that the mean energy of H' for the states \1'' of Ti' coincides 
with the mean energy of H for those states of Ti obtained through the embedding T, namely, 

< '^'\H'\m' >=< T^'\H\T^' > (8) 

In other words T\I'' is used as a variational state for the eigenstates of the Hamiltonian 
H. In particular T should be chosen in such a way that the states truncated in 7i , which 
go down to H', are the ones expected to contribute the most to the ground state of H. 
Thus Eq. (^ is the basis of the so called variational renormalization group method (VRG) 
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As a matter of fact, the VRG method was the first one to be proposed. The ERG came 
afterwards as a perturbative extension of the former (see later on). 

More generally, any operator O acting in Ti. can be "pushed down" or renormahzed to a new 
operator O' which acts in Ti' defined by the formula, 

O' = T^OT (9) 

Notice that Eq.(|^) is a particular case of this equation if choose O to be the Hamiltonian 
H. 

In so far we have not made use of the all important concept of the block, but a practical 
implementation of the VRG or ERG methods does require it. The central role played by 
this concept makes all the real-space RG-methods to be block methods. 

Once we have established the main features of the RG-program, there is quite a freedom 
to implement specifically these fundamentals. We may classify this freedom in two aspects: 

• The choice to how to reduce the size of the lattice. 

• The choice of how many states to be retained in the truncation procedure. 

We shall address the first aspect now. There are mainly two procedures to reduce the size 
of the lattice: 

• by dividing the lattice into blocks with ns sites each. This is the blocking method 
introduced by Kadanoff to treat spin lattice systems. See figure 1. 

• by retrieving site by site of the lattice at each step of the RG-program. This is the 
procedure used by Wilson in his RG-treatment of the Kondo problem. This method 
is clearly more suitable when the lattice is one-dimensional. 

We shall be dealing with the Kadanoff block methods mainly because they are well suited 
to perform analytical computations and because they are conceptually easy to be extended 
to higher dimensions. On the contrary, the DMRG method introduced by White works 
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with the Wilsonian numerical RG-procedure what makes it intrinsically one-dimensional and 
difficult to be generalized to more dimensions. This situation has changed recently in part 
as S.R. White has devised a numerical improvement of the DMRG which is applicable to a 
1/5-depleted 2D lattice We have formulated our Variational and Fokker-Planck DMRG 
procedures as block renormalization methods [O. 

Block RG-methods have recently received also renewed attention in one- dimensional prob- 
lems in connection to what is called a quantum group symmetry [T^, |jl6|. Based upon 
this symmetry we have constructed a new BRG-method that we call g-RG which among 
other features it is able to predict the exact line of critical XXZ models in the Anisotropic 
Heisenberg model, unlike the standard BRG-method. 



III. THE ROLE OF BOUNDARY CONDITIONS AND REAL-SPACE RG 

The ffist advance in trying to understand the sometimes bad numerical performance of 
the BRG methods came in the understanding of the effect of boundary conditions (BC) on 
the standard RG procedure p^ . 

White and Noack [|l^ pointed out that the standard BRG approach of neglecting all 
connections to the neighbouring blocks during the diagonalization of the block Hamiltonian 
Hb introduces large errors which cannot be corrected by any reasonable increase in the 
number of states kept. Moreover, in order to isolate the origin of this problem they study an 
extremely simple model: a free particle in a ID lattice. As a matter of fact, it was Wilson 
|[T8| who pointed out the importance of understanding real-space RG in the context of this 
simple tight-binding model where the standard BRG clearly fails as we are going to show. 
The reason for this failure can be traced back to the importance of the boundary conditions in 
diagonalizing the states of a given block Hamiltonian Hb in which the lattice is decomposed 
into. Notice that in this fashion we are isolating a given block from the rest of the lattice 
and this applies a particular BC to the block. However, the block is not truly isolated! 
A statement which is the more relevant the more strongly correlated is the system under 
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consideration. Thus, if the rest of the lattice were there it would apply different BC's to the 
boundaries of the block. This in turn makes the standard block-diagonalization conceptually 
not faithfully suited to account for the interaction with the rest of the lattice. 

Once the origin of the problem is brought about the solution is also apparent: devise a 
method to change the boundary conditions in the block in order to mimick the interaction 
with the rest of the lattice. This is called the Combination of Boundary Conditions (CBC) 
method which yields very good numerical results. This method has not yet been generalized 
to interacting systems. However in reference an alternative approach is proposed under 
the name of Density Matrix Renormalization Group (DMRG) which applies to more general 
situations and also produces quite accurate results. 

In a recent paper |21]] we have reconsidered again the role of BC's in the real space RG 
method for the case of a single-particle problem in a box. The continuum version of this 
Hamiltonian is simply H = —-§^- We shall consider open chains with two types of BC's at 
the ends: 



Fixed BC's: ipiO) = ^{L) = 



(10) 



Free BC's: ^{0) = ^{L) = 
ox ox 



(11) 



The lattice version of H for each type of BC's is given as follows: 



/ 



Fixed 



2 -1 
-1 2 -1 
-1 2 



V 



2 -1 
-1 2 



/ 



Hpree 



1 -1 
-1 2 -1 
-1 2 



V 



2 -1 
-1 1 



(12) 



The only difference between Hpixed and Hpree appear at the first and last diagonal entry 
(2 ^ 1). The exact solution of (|12D is very well-known and we give it for completeness: 
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Fixed EC's: Mj) = ^^^/ j, = 4sin2(^A__^) (13) 



7r?7 1 TTT? 

Free EC's: ^.(j) = ^ cos — (j - -), i^^. = 4 sin^l— ) (14) 



J = l,2,...,iV; n = 0,l,...,N-l. 

where the N^^s are normahzation constants and is the number of sites of the chain. 

Eefore getting into the problem of the renormahzation of these Hamiltonians, it is worth 
to pointing out another physical realization of Hpree'- A simple magnon above a ferromag- 



netic background satisfies Free BC's. See I^Tl for more details. 

Now let us get to the problem of renormalizing the tight-binding Hamiltonians (0). 

In figure 2 we show the ground state and first excited states of the chain with fixed and 
free EC's. 

It is clear from Fig. 2 that a standard Block RG method is not appropiate to study the 
ground state of fixed BC's since this state is non-homogeneous while the block truncation 
does not take into account this fact. 

Each piece of the ground state within each block satisfies EC's which vary from block to 
block. This is the motivation of reference |]T^ to consider different EC's in the block method, 



yielding quite accurate results. 

We observe that the standard RG method performs rather poorly as compared to the CEC 
method which yields quite the exact results. 

The other alternative to the CEC method is the Density Matrix RG method which can be 
phrased by saying that the rest of the chain produces on every block the appropiate EC's 
to be applied to its ends, and it has the virtue that can be generalized to other models, 
something which is not the for the CEC method. 

On the other hand, the ground state of Hpree is an homogeneous state (see Fig. 2) which 
in turn suggests that a standard RG analysis may work for this type of EC's. We shall 
show that this is indeed the case if the RG procedure is properly defined. The key of our 
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RG-prescription is to notice that Hpree has a geometrical meaning: Hpree is the incidence 
matrix of a graph, and it is called minus the discrete laplacian —A of that graph. Notice 
that Hpixed has not such geometrical interpretation, in fact, it concides with the Dynkin 
diagram of the algebra A^. 

Based on this observation the Kadanoff blocking is nothing but the breaking of the graph into 
N/us disconnected graphs of sites each. We shall choose = 3 in our later computation. 
The previous geometrical interpretation of Hpree suggests that we choose the block 



Hamiltonian Hb to be the incidence matrix of a disconnected graph, namely, 

/ 



1 -1 
-1 2 -1 
-1 1 



1 -1 
-1 2 -1 
-1 1 



V 



■ / 











1 -1 
-1 1 



1 -1 
-1 1 



V 



(15) 



andthe interblock Hamiltonian Hbb above describes the interaction between blocks. 
Hbb in turn also coincides with the incidence matrix of a graph which contains the missing 
links which connects consecutive blocks. In a few words: our RG-prescription introduces 
free EC's at the ends of every block. This condition fixes uniquely the breaking of Hpree 
into the sum Hb + Hbb- This is the choice we make. It should be emphasized that the 
splitting of Hpree i^to two parts Hb + Hbb is by no means unique, so that different choices 
may lead to very different results. 

Prior to any computation we notice that the previous RG-prescription should lead to an 
exact value of the ground state energy, for the ground state of each block is again a constant 
function. The question is therefore to what extent our method is capable of describing the 
excited states. We shall concentrate ourselves to the first excited state since computations 
can be carried out analytically. 
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First of all we diagonalize Hb within each block of 3 sites, keeping only the ground 
state ■0o'^^ and the first excited state ipf'^ (3^2 truncation). The superscript denotes the 
initial step in the truncation method. In figure 3 we picture the 3 eigenvectors for the 3-site 
Hamiltonian which will be the building blocks for our BRG, namely the two lowest ones. In 
the standard RG method we would choose V'o"^ and if^f^ as the orthonormal basis for the 
truncated Hilbert space and obtain the effective Hamiltonian H'^ and H'^^. In our case it 
is convenient to express these effective Hamiltonians in a basis expanded by the following 
linear combination: 



(0) 



1 



(0)s 



(16) 



which are also an orthonormal basis of the truncated Hilbert space. 
In this basis the truncation of Hb reads as follows, 



H'b- 



A 
0^0 
A 



\ 



A 



A 



1 1 -1^ 



-1 1 



(17) 



(18) 



V / 

with e taking on the value e*^"-* = 1 in the initial step of the RG-method, which is the energy 
of the state ■ 

The truncation of Hbb is more complicated, the result being: 
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H 



BB 



BB 



B C 
C D + B C 
C* D + B 



D + B C 
C* D 



(19) 





( 2 A 

a ab 




' -ab -a" \ 




ab^ 


B = 


^ab b^ ^ 


C = 


^ —b"^ —ab ^ 


D = 


^ab a"^ J 



with a b taking on the values a*-"-* = ^ ~ ^ aiid b^^^ = ^ + | in the initial step of the 
RG-method. 

The nice feature about the basis (p!^)-([T7|) is that all rows and columns of (|1^) and (pJ]) 
add up to zero, just like the original Hamiltonians (|1^), implying that the constant vector 
is an eigenvector with zero eigenvalue of the renormalized Hamiltonian! 

We shall call HN/^{e, a, b) the sum of the Hamiltonians (0) and (|19]) for generic values of 
e, a and b. Next step in our RG-procedure is to form blocks of 4 states of the new Hamiltonian 
His[/z{e, a, b) and truncating to the two lowest ip^^ and ip^i^ energy states within each 4-block 
(4 — i> 2 truncation). The reason for this change in the number of sites per block (from 3 to 



4) is motivated by the form of H'^^ in (|T9D and the fact that if we try to make a second 
step in the RG-method with 3-blocks the method is doomed to failure because the constant 
state of Fig. 2 would no longer be the ground state. 

Fortunately enough, with 4-blocks if we define new states and ip^l'^ in the same form as 
we did in Eg. (p!2D , we obtain that the new effective Hamiltonian is obtained by a redefinition 
of the parameters, namely. 



HN/3{e, a, b) — ^ Hjy/G^e', a', b' 



- + a' + b'-A 



(20) 
(21) 
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2^2 



a + 6 — 



a[a 



36^ + A) + f 
A(A + a2 - 62) 



(22) 



6' 



2^2 



a + b 



a{a^ - 362 + A) + I 
A(A + a2 - 62) 



(23) 



A 



- 62)2 + (_ _ 2a6)2 



(24) 



In this fashion, the constant state of Fig. 2 is again the ground state of the model and 
moreover, upon iteration of Eqs.(|20|)-(p^ there are no level crossing among the excited 
states. Otherwise stated this means that the level structure of the block Hamiltonian Hb is 
preserved under the action of our BRG-method based upon the reduction from 4 to 2 states. 

The energy Ei{N) of the first excited state of a chain with = 3 x 2"^ sites can be 
obtained iterating m times Eqs.(|20| ) - (|2^) : 

(25) 



Ei{N = 3x 2™) = e^'") 



The initial data are given by: 

^(0) 



76 2' V6 2 



(26) 



For low values of N the deviation of e^"^) with respect to the exact result is small (see [pl[|). 
Recall that we are only keeping two states in our RG-procedure, and that the ground state 
energy is exactly zero by construction!. But what is more interesting about these (see pT| ) 
is that we are able to obtain the correct size dependence, i.e., l/N"^ of e^"^\ As a matter of 
fact, the energy of the first excited state behaves for large N as ([I^: 



Et^^'\N) c,^^jN\ with c, 



exact 



(27) 



while our BRG-method gives. 



Ei^'^^V) ~ cbrg/N\ with cbrg = 12.6751 



(28) 
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The achievement of the is a remarkable result which in turn allows us to match 

the correct order of magnitude of the energy. For instance, for 10 iterations our RG-method 
with 2 states kept gives the energy of the order of 10~^, which is precisely the same order of 
magnitude as for the CBC method but with 8 states kept in the case of Fixed BC's. Recall 
that the standard BRG performs as bad as a 10~^ order of magnitude. 

IV. WAVE-FUNCTION RECONSTRUCTION 

Insofar we have only used our RG method to "reconstruct" the energies of the lowest 
states step by step, but we can also use this method to reconstruct the shape of the associate 
wave function in the real space. This simple fact leads to the consideration of RG applications 
beyond the original scope for which it was devised. As in this fashion we are making a picture 
of the wave function, it is natural to use the RG as a image compression method for coding 
images in order to facilitate its transport through the networks. For more details see the 
notes by Nishino in these proceedings |Q. 

We are able to make a reasonable picture of the first excited state wave-function based 
upon our BRG-procedure when compared with the exact form depicted in Fig. 2. As we are 
working with a real-space realization of the renormalization group method, this is something 
we have at hand. To do this we need to perform a "reconstruction" of the wave-function. 
This reconstruction amounts to plot the form of our aproximate wave-function in each and 
every of the 3-blocks out of the 2"^+^ in which the original chain is decomposed into under 
the BRG-procedure. Recall that in the initial step we started out with blocks of 3 states 
keeping the two lowest states i/jq'^ and ipi^^ (3 — 2 truncation). In the next step we make 
blocks of 4 states keeping the two lowest states ipl)'^ and (4^2 truncation) and then we 
perfom the iteration procedure over and over. As a result of this procedure we may express 
the two lowest wave functions of the m + 1-th step in terms of those of the previous m-th 
step by means of the following matricial form: 
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(m+l) 



V2 



J 



(m) 



1 

72 



V 



I 



(29) 



L \ / \ / R 

where the LHS of Eq. (^) represents the wave function of 3 x 2"^+^ sites while in the RHS 
we have a left-wave-function of 3 x 2™ sites and another right-wave-function of 3 x 2™ sites, 
so that everything squares. The parameters appearing in Eq. (E9|) turn out to be given by: 



(30) 



with Am as in Eq, 
We may recast Eq. 



^ — On h 1 



(al - bl + A„) 



2VAm(Am + a^-6L 



Pn 



Their initial values are ao = 1/vlO and Po = 3/VlO. 
) in more compact form by writing: 



(31) 



^(m+l) ^ ^ ^(m) ^ ^(m) 



(32) 



where 



1 ( 1 



72 



R„ 



J 



1 

V2 



1 



(33) 



We may call Eq. (p^) the reconstruction equation. This is the master equation that when 
iterated "downwards" (reconstruction) allows us to obtain the picture of our approximate 
BRG-wave-function corresponding to every and each block of 3 sites of the 2™'^^ blocks in 
which the chain is decomposed into. At the end of the iteration procedure we end up with 
expressions for the values of the 3-sites wave-functions in terms of the initial two lowest 
states ipl)'^ and The first one is a constant function while the second is a straight 

line of negative slope. Thus, these two states turn out to be the building blocks of our 
B RG-pro cedure . 

When using the reconstruction equation to obtain the wave function we may use a 
binary code based upon the labels L (left) and R (right) to keep track of the different 3-sites 
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blocks which make up the chain. Thus, in one dimension the RG-blocks are in a one-to-one 
correspondence with a binary numerical system. In general, for other dimensions we may 
state squematically the following correspondence: 

BRG-prescription < — > "Number System" 

This simple observation is the basis of a coding system for compressing pictures whatever 
may be its origin. In two dimensions we need more digits to make the coding, but it works 
likewise and serves for image compression @]. 

V. THE CORRELATED BLOCK RENORMALIZATION GROUP (CBRG) 

We have already mentioned in the Introduction that the DMRG is not only a powerful 
computational method but also a source of inspiration for further works concerning the RG. 
For these reasons, it is worthwhile to explore different options or alternatives to the DMRG 
which may be useful in situations where the DMRG encounters difficulties, as in the case of 
2D quantum systems. The main message of the DMRG is that blocks are correlated. The 
implementation of this idea by means of the density matrix formalism may be not the unique 
way to proceed. On the other hand, the "onion-scheme" a la Wilson adopted by the DMRG, 
while being one of the reasons of its spectacular accuracy, imposes certain limitations. 
At this stage it is not clear how fundamental are the density-matrix formalism or the onion- 
scheme for a RG method which takes into account the correlation between blocks. One can 
indeed combine the Kadanoff block method with the use of a density matrix in the process 



of truncation, as in reference |12|. More work remains to be done to see wheather there is a 



real improvement of the standard BRG method by combining it with the DMRG as in [jT2[ 



In this section we want to explore another possibility which is to give up both the density 



matrix and the onion-scheme (see [^). With this point of view in mind, it would seem that 
we should come pretty close to the standard BRG method, were it not for the enormous 
freedom hidden in a Real-Space RG method. This freedom comes from the separation of the 
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Hamiltonian into an intrablock Hb and an interblock Hbb Hamiltonian. This is a source of 
ambiguities which can be sometimes mitigated with the aide of symmetry arguments, but 
not fully eliminated though. This ambiguity shows up specially for terms in the Hamiltonian 
acting at the boundaries of the block. There are no general criteria as to how to include 
this type of terms either into the intrablock or into the interblock Hamiltonians, or into 
both! For example, in the ID Ising model in a transverse field (ITF model), a choice which 
preserves the selfduality of the model attributes some self-couplings to the Hb and others 
to the interblock Hbb, and it yields to an exact value of the critical point and the critical 



exponent z/ [^, ||T5[. The ambiguity in the splitting of H into the sum Hb + Hbb thus affect 
deeply the truncation procedure itself, which is based on the diagonalization of Hb- Rather 
than blaming the BRG for its lack of uniqueness, we should use its freedom to allow the 
blocks to become correlated in the RG procedure. In our present approach this correlation 
will be taken into account in a "dynamical" way rather than in a "statistical" way as in the 
DMRG. This will be achieved by the introduction of interblock operators which reflect the 
"influence" between neighbour blocks and which are defined at the boundary of the block 
in the first step of our CBRG method. 

We have chosen to illustrate our approach the ID and 2D tight-binding models mainly 
for simplicity reasons, but we believe that our method could be applied to more complicated 
problems. In fact, the first step in this direction was already undertaken in reference ||21|| , 



where only 2 states at each stage of the RG-blocking were retained. This in turn allowed us 
to obtain the l/N"^ scaling law for the size dependence of the first-excited-state energy. 



We shall give the general mathematical structure underlying the results of reference [21 



This will allow us to retain more than two states in the RG-truncation and also to consider 
the two-dimensional tight-binding model. In this fashion, we shall recover the n'^/N'^ scaling 
law for the n-th excited state of the ID model and the scaling law ""^^2^ in the 2D case. 
These results will then show that the CBRG method describes correctly the low energy 
behaviour of the ID and 2D Laplacian. 
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VI. THE CBRG METHOD: ONE DIMENSION 



The problem we want to study is the one-dimensional Tight-Binding model in an open 
chain with different boundary conditions at its ends. The Hamiltonian for this system takes 
the following matricial form, 



Hb,b' 



b -1 
-1 2 -1 
-1 2 



V 



2 -1 
-1 b' 



(34) 



where b and b' take on the values 1 (or 2) corresponding to Free (or Fixed) BC's respectively. 
This Hamiltonian is the discrete version of the Laplacian H = —d"^, while the Free or Fixed 
BC's correspond in the continuum to the vanishing of the wave function (Fixed BC's) or its 
spatial derivative (Free BC's) at the ends of the chain, i.e.. 



b = 2 ^ ^(0) = Fixed BC 
6 = 1^ 11(0) = Free BC 



(35) 



and similarly for b' which contains the BC at the other end of the chain. 
Hence, altogether there are 4 Hamiltonians of the type in (|5^, whose eigenstates and eigen- 
values are the subject of our RG-techniques. 

The first step in the RG method is to divide the lattice into blocks containing sites 
each and labeled with and index p (= 1, . . . , N/ug). Let us suppose for a moment that we 
isolate the pth-block from the rest of the lattice so that its dynamics, as an independent 
entity, is governed by a Hamiltonian denoted by Ap, which we may call uncorrelated block 
Hamiltonian. The restoration of the block back into the lattice involves two effects. The 
first one is that the BC's of the p-th block may change under the influence of the p + 1 and 
p — 1 blocks. We describe this change of BC's by the action of Boundary Operators denoted 
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by Bp^p±i on the pth-block. The second effect is the interaction between the pth-block and 
its neighbours p + 1 and p — 1, given by interaction Hamiltonians Cp,p±i which act on both 
p and p + 1 blocks simultaneously. If the problem under consideration is translationally 
invariant, all the Hamiltonians defined above are independent of the block label p, in which 
case we denote them by, 

Ap = A 

Bp,p+i = Br -Bp,p_i = Bl (36) 
Cp,p+i = C Cp^p^i = 

The Hpree,FTee Hamiltoniau (p^ gives an example of this as we shall show below. Hence, 
for the time being, we shall consider the situation described by (|36| ) and leave the more 
general case after explaining the general ideas. 

In the standard BRG method the block Hamiltonian Hb and the interblock Hamiltonian 
Hbb are given, according to our previous definitions, by the following formulas 

Hb = A + Bl + Br (37) 



H 



BB 



(38) 



The whole Hamiltonian is by all means the sum of Hb and Hbb for all the blocks of the 
chain. 

Next step in the RG method is to diagonalize Hb and keep its, say m {m < Ug), lowest 
eigenstates. The truncation is given by a x m matrix T whose columns are precisely the 
components of the m lowest eigenstates of Hb- The renormalized Hamiltonian in the new 
basis is given by, 

H' = T\Hb + Hbb)T (39) 
At first sight from Eq. (|57D it would seem that we have taken into account the effect of 



the EC's on a given block. However, as the examples show, this is quite a bit illusory. 
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On the other hand, the distinction among A, Bl and Bn is rather inmaterial as far as 
Hb is concerned, and in fact no distinction of this sort is made in the standard BRG 
formahsm. Finally, let us observe that Hb and Hbb play rather different roles in the 
truncation procedure. This asymmetry has been observed as a source of problems by several 
authors in the past 



We shall mention that this asymmetry has recently been related to quantum groups in 
a fashion which has led to a new RG method called the Renormalization Quantum Group 
method [|l5l, ||T6 |. 

Therefore, from various points of view, one is urged to make more explicit the role 
played by the BC-operators B^ and Bji in our CBRG procedure. For this purpose, we have 



found convenient to use the concept of superblock already introduced in reference [T^. We 
shall define a superblock as the set of two consecutive blocks, p and p + 1 and denoted by 
{p,p + 1). The great advantage of the superblock is that it allows us to materialize the 
distinction among A, Bl and B^ . In fact, just as the isolation of a single block leads us to 
the definition of the Hamiltonian A, the isolation of two blocks contained in a superblock 
allows us to define B^, Bji and also C through the superblock Hamiltonian HgB as follows. 



B 



V 



(40) 



A + Br c 

A + Bl 

Similarly, the Hamiltonian describing the interaction between superblocks is given by (see 
Fig.4) 



sB,sB 







Br C 

Bl 



\ 







(41) 



Now instead of diagonalizing Hb in Eq. (|37D, in the CBRG method we shall diagonalize H^b 
in Eq. (^OD, and afterwards keep the m = Ug lowest eigenstates in the tight-binding model. 
As in the standard BRG method, the change to the truncated basis defines the renormalized 
operators as follows: 
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sB 



T^H.rT = A' 



(42) 



H. 



sB,sB 



B,sB 



T 



C'^ B' 



(43) 



where the matrices A', 5^, B'j^ and C are the renormahzed version of the operators A,Bji, 
Bl and C, and they exhibit the same geometrical interpretation for the renormahzed block 
as their unprimed partners for the original blocks. 

If we set Br = Bl = in Eqs. (|iD| ) and (0), then after the first RG-step we get 
5^ = B'l^ = and thus the previous RG-scheme coincides with the standard BRG. We may 
say that uncorrelated blocks are in a sense a fixed point of our method. However, this fixed 
point may be unstable, and to explore this possibility one has to look for non-vanishing 
5-operators and their RG-evolution. 

Let us address now some examples. We shall first study the Hamiltonian (|3^) with Free 
BC's at the ends {b = b' = 1). Choosing = 3 for example, we see that the choice for the 
operators A, Br, Bl and C in the first step of the CBRG procedure is given by. 



A 



1 -1 
-1 2 -1 
-1 1 



Bf 







V 



1 



/ 



B, 



\ 



V 







0^ 



c 



V 




-10 



(44) 



This choice is equivalent to the assumption that an isolated block satisfies Free BC's at its 
ends. The role of Br and Bl is to join these blocks into a single chain. This is the geometrical 
explanation of Eqs. (0). In more general cases one must have to explore which is the best 
choice. The generalization of Eqs. (|4^) to blocks with more than 3 sites is obvious. In Table 
1 we collect our CBRG-results for the first 5 excited states for a chain of = 12 x 2^ = 768 
sites. Comparison with exact results gives a good agreement. 

An important feature of our CBRG method is that the n'^/N'^- scaling law (A^ — > oo) 
for the energy of the n-excited states of a chain made up of A^ sites, is reproduced correctly 
(see Fig. 5). In Table 2 we show the variation of the first-excited-state energy with the size 
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of the chain. From those values we can extract the corresponding 1/A^^-law which turns 
out to be, 

Ep««)(iV) = cg^G^, cgi^G = 9.8080, (N oo) Free-Free EC's (45) 

while the exact value for the proportionality constant c is Cg^act = = 9.86. This amounts 
to a 0.6 % error. 

Likewise, we have enough data so as to obtain the corresponding n^/A^^-law for the whole 
set of 5 excited states. Thus, the scaling law we obtain is, 

2 

^«^)(iV) = ccBRG^, ccBRG = 8.4733, (N oo) Free-Free EC's (46) 

which now amounts to a 7.34 % error. This is a natural fact from the worse knowledge of 
the highest excited states of the spectrum in a RG-scheme. 

We can make even more explicit the successful achievement of the 1/A^^-scaling law by 
leaving as a free adjustable parameter the exponent of in addition to the proportionality 
constant. Let us denote by 9 this critical exponent. Using data from 20 to 50 steps of our 
CERG-method for several truncation of states according to our scheme 2ns (namely, 
Us = 10, 13, 20) we arrive at the following results, 

e(cbrg)^j^^ = ccBRG^, {N oo) Free-Free EC's (47) 
For 20 — ^ 10 9 = 1.9708 

For 26 — ^13 9 = 1.9734 (48) 
For 40 — ^20 9 = 1.9854 

These results clearly support the fact that we have correctly reproduced the exact value of 
9 = 2 for the finite-size critical exponent. 

Last, but not least, as was proved in |^T| our CERG method gives the exact energy of the 
ground state for every step of the RG-procedure for Free-Free EC's. 

In tables 1 and 2 we also show the results we have obtained with a DMRG analysis following 
White's method DTj. This analysis is based on the onion scheme of enlarging the lattice site 
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by site a la Wilson. The results coincide with the exact values within the 4 digits precision 
used here, but they start differing when keeping more digits. Nevertheless, the DMRG is 
much more time consuming than our CBRG method for it has to build the lattice site by 
site, while the CBRG reproduces the lattice by blocking which is much more efficient as far 
as CPU time is concern, and moreover, it applies to two-dimensional situations where the 
onion scheme fails to reproduce the lattice. We have also performed the DMRG analysis in 
ID for Fixed-Fixed EC's in table 4 where the same considerations apply. The CBRG method 
is also more suitable for analytic formulations [^. 

In reference it was shown that one can reproduce easily the wave function of the excited 



states. This procedure was called reconstruction since it works "downwards" in the CBRG 
method. The basic equation to be used is the reconstruction equation pi] . 



^(^•+1) = lMI^ + rM^ (49) 

where denotes the collection of m lowest eigenstates in the r-step of the CBRG- 
procedure, and L,,, Rr are the block matrices in terms of which the truncation matrix 

can be written as T'^ = {Lr, Rr)- 
Our results for a chain of = 12 x 2^ = 768 sites and = 6 states kept are given in Fig. 6 
where we have plotted the first 5 excited states and compare them with the exact wave 
functions. There are some remarkable facts regarding these figures. Firstly, the number of 
nodes is correctly preserved by our CBRG wave functions. Secondly, the Free-Free type of 
boundary conditions are also correctly reproduced at the ends of the chain. And lastly, it is 
worthwhile to point out that the CBRG wave functions "degrade gracefully" as the energy 
of the excited state raises in accordance with the fact that the lower the energy is, the more 
reliable are the results. 

This ends the results for the Free-Free BC's. In order to address other types of BC's we 
must come back to the case where the matrices A,B{i, Bl and C depend on each particular 
block. Thus, for example, for the Fixed-Free BC's we shall choose as the uncorrelated 
A-matrix for the block located to the left end of the chain the following form [ug = 3), 
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^2 -1 ^ 



^1 



V 



-1 2 -1 
-1 1 



Fixed-Free EC's 



(50) 



while the remaining matrices Ap, (p = 2, . . . , A^/3), will be given by Eqs.(P^), (0). 

For Free-Fixed EC's, it is the last A-matrix which we have to take different from the 
others, namely, 



/ 



A 



1 -1 
-1 2 -1 
0-12 



\ 



Free-Fixed EC's 



(51) 



As for the Fixed-Fixed EC's case, we must change the A-matrix at both ends of the chain 
according to the following prescription. 



/ 



V 



2-10 
-1 2 -1 
-1 1 



A 



N/3 



V 



1 -1 
-1 2 -1 
0-12 



Fixed-Fixed EC's 



(52) 



Then we follow the same steps as for the Free-Free EC's, taking care that the A,Bji, and 
C matrices in each CERG-step may depend on the position of the blocks. This implies in 
particular that the embedding T-matrices may also vary from block to block. 

In Tables 3 and 4 we summarize our results for the Free-Fixed and Fixed-Fixed EC's 
(Fixed- Free EC's are equivalent to Free-Fixed EC's by parity transformation). In these 
tables we present our CERG results for the first 6 lowest lying states for the ID tight- 
binding model in a chain of = 12 x 2^ = 384 sites with mixed boundary conditions, 
and they are compared against the exact and standard ERG values. Several remarks are in 
order. First, we observe that the CERG method produces a good agreement with the exact 
results and certainly much more accurate by several orders of magnitude than the old ERG 
method. Second, the CERG method is able to reproduce the corresponding n^/A^^-scaling 
laws for the spectrum of excited states in each case of mixed EC's. Namely, 
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• For Free-Fixed BC's and considering just the ground state, we have 

EiCBRG)^^^ = cP^nc^. cgLc = 9.072, (N oo) Free-Fixed BC's (53) 

which amounts to a 8 % error with respect to the exact value of Cexact = tt^. 
As for the corresponding law for the whole spectrum, we find 

^^^)(iV) = ccBRG ^'^^^l^\ ccBRG = 7.6729, (AT ^ oo) Free-Fixed BC's (54) 

which represents a 11.5 % error with respect to the exact value of vr^. 

• For Free-Fixed BC's and considering just the ground state, we have 

^^^)(iV) = cS^G^, cP^RG = 8.35, {N oo) Fixed-Fixed BC's (55) 

which amounts to a 8 % error with respect to the exact value of Cexact = tt^. 
As for the corresponding law for the whole spectrum, we find 

EiOBRG)^^-^ = ccBRG ^''^P\ ccBRG = 6.9696, (A^ oo) Fixed-Fixed BC's (56) 

which represents a 16 % error with respect to the exact value of vr^. 

We obtain bigger errors in the determination of these scaling laws as compared with the Free- 
Free case mainly because we have used less data in our fitting. Nevertherless, we find a good 
agreement with the exact results. Yet, there is another reason as to why the accuracy in the 
case of mixed BC's is worse, namely, the ground state wave function \E'o is not homogeneous 
in space as it is in the Free- Free case |2T|. This makes the RG-procedure more involved and 



a source of extra uncertainties. 

Let us mention in passing that we are also able to make a wave function reconstruction in 
the mixed BC's cases as has been done for the Free-Free BC case. 

The outcome of all the results presented so far is that we have succeded in devising a 
Real-Space RG method capable of reproducing the correct eigenvalues and eigenstates for 
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the tight-binding model as originally envisaged by Wilson, within a certain accuracy which 
can in principle be improved. 

Althoug the model we have employed to test our CBRG-method is a tight-binding model, 
there are some remarkable facts regarding the fixed-point structure of our CBRG-solution 
that we would hke to stress. Namely, we have found that after enough number of CBRG- 
iterations, the matrices A, Bl, Br and C in the Free- Free case scale nicely with the size 
N of the chain according to the dynamical critical exponent z. To be more precise, let us 
introduce the fixed point values of those matrices denoted by A*, Bl, i?^ and C* which we 
define as, 

A* = N-'a*, Bl = N''bl, Bl = N-'b% C* = N-'c*, Fixed-Point values (57) 

in terms of the scaled matrices a*, bl, bl and c*. For a block of 3 sites (n^ = 3) we find 
the following Fixed- Point structure parametrized by two constants s and t (for bigger we 
need extra parameters). 



(58) 
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with s — 1.3993 and t — 1.9581. The critical exponent z we obtain is, 

z = 0.9999 



(59) 



which is indeed very close to the exact value z = 1 (actually, it differs in the ninth decimal 
digit). 

The interpretation of this Fixed-Point in the context of the CBRG method is as follows. We 
pointed out before that when the boundary operators vanish we recover the standard 
BRG method in which the blocks are not correlated. Here we find that it is the uncorrelated 
Hamiltonian which vanish, while the boundary 6^ r and interaction c operators do not vanish 
within the scaling law. This fact may perhaps be interpreted by saying that in the example 
under study the correlation between blocks is more important than their selfenergy. In 
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references , , ||2^ it was shown that the DMRG method leads, in the thermodynamic 
hmit, to a "product form" ansatz for the ground state wave function. In our case we see from 
Eqs.([57|), (|58| ) that we also reach thermodynamical limit, which leads us to ask about the 
nature of the ansatz for the ground state and excited states implied by the CBRG method. 
The answer to this question will be addressed in a future publication but it suffices to say 
that both the DMRG and the CBRG methods seem to yield different ansatzs of the ground 
state wave function. In a few words, the DMRG is associated with a "vertex picture" while 
the CBRG is associated with a "string picture" . 



VII. THE TWO-DIMENSIONAL CBRG- ALGORITHM 

The RG-method that we have devised in the one- dimensional problem can be generalized 
in a natural way to higher dimensions. We shall consider for simplicity the 2D case. First 
of all, we divide the square lattice into blocks of sites each. Each block will in turn be 
a square lattice with a minimum of 4 sites (=2x2 block). As in ID, we shall define the 
following Hamiltonians to carry out the CBRG-program, 

• Ap = self-energy of the p-th block isolated from the lattice. 

• Bp^q = self-energy of the p-th block induced by the presence of the g-th block. 

• Cp^q = interaction between the p-th block and the g-th block. 

The difference with respect to the ID case is that each block has now 4 neighbours and 
therefore there are four different B and C matrices. 

Let us consider again the Hamiltonian of a free particle moving in a 2D-box with Free 
BC's at the boundaries of the box. The 2D Hamiltonian is given again by the incidence 
matrix of the lattice. As in ID we shall choose the matrix A as the incidence matrix of the 
block. Thus, for example, for a 2 x 2 block we have. 
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A 



^2 -1 -1 ^ 



-12-10 
0-12-1 
-10-12 



(60) 



The 4 Boundary Operators B have a diagrammatic representation which helps us to 
keep track of their location in the block Hb and interblock Hbb Hamiltonians. Their explicit 
matricial form is as follows, 



Bi2 = B43 = Bl 







B21 = B34 = Br 



\ 



(61) 



B 



14 



B 



23 



B 



D 



\ 



1 



B, 



41 



5. 



32 



B, 



u 



\ 







(62) 



where the labels denote the position of the neighbouring blocks and we have used the trans- 
lation invariance of the 2D tight-binding model so that we need only to distinguish between 
Right and Left, and Up and Down. 



As for the Interaction C-Operators ||2^ we have the following matricial representation, with 
the same considerations as for the -B-operators, 



C12 — C. 



43 



c 



LR 



^ ^ 



-10 
-1 




C21 — C:- 



34 



c 



RL 



^0-1 0^ 





0-10 



(63) 



29 



C^A — C; 



23 



C 



DU 



^ 0^ 




0-100 
-10 



C41 — C: 



32 



UD 



-1 
0-10 





(64) 



Thus translation invariance reduces the number of independent CBRG-matrices by a half. 
These relations are particular of the problem at hand but we must left open the posibility 
of having all those matrices different from each other in order to handle more complicated 
problems. 

Now that we have all the elements entering in our CBRG-method we proceed to construct 
the block HgB and interblock if^s sB Hamiltonians out of them. To this end we have to 



consider a superblock made up of 4 blocks |22[. Thus, for HgB we have 



B 



A + Bl + Bd Clr 
Cbl a + Br + Bd Cdu 

Cud A + Br + Bu 



Cdu 


Crl 



\ 



c 



UD 







c 



LR 



(65) 



A + Bl + Bu 

This is a Aug x matrix made up of x rig matrices. 

As for the interblock Hamiltonian Hsb,sB we have to distinguish between (s-B, S-B)-couplings 
of horizontal type denoted by H^^^^^ and vertical type denoted by H^^b^Jbi which read 
explicitly as, 



30 



H. 



{hor) 



sB,sB 



Br Crl 





Bn Crl 


ClrO Bl 

Clr Bl 





(ver) 
sB,sB 



B 



D 



Ci 



Bd Cud 







Cz5[7 Bu 
Cdu 



UD 







Bu j 
(66) 



where we have made use again of translational invariance. 

Once that we have made our choice for the decomposition of the total Hamiltonian of the 
2D-tight-binding model into block and interblock Hamiltonians according to our CBRG- 
prescription, we can carry on with the truncation part of the RG-method. We shall keep Us 
states out of states per superblock so that our truncation scheme may be summarized 
as, 

4ns (superblock) — > Ug (new block) 

Recall that at each step of the CBRG-method we need to identify the A, B^, Br and 
C operators which define the truncation procedure for the next step of the method. For 
this purpose, firstly the truncation of the superblock HgR gives rise to the A' uncorrelated 
self-energy operator for the next RG-step, namely. 



HsB ( 4ris X 4ns matrix) — > A' (n^ x Ug matrix) 



(67) 



To identify the rest of the operators we have to renormalize the interblock Hamiltonian which 
comes in two types, horizontal and vertical. The renormalization of the i^s^^Is Hamiltonian 
is given by p2| . 



H. 



(hor) 



sB,sB 



^ B'r C'^l"^ 
\C'lr B'^ 
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(68) 



Likewise, for the H^sb^]^ Hamiltonian we have, 



H 



(ver) 
B,sB 



B'd C'ud ^ 



y C'du J 



(69) 



Now that we have identified all the operators defining the CBRG method at the new stage 
of the renormalization, we may reconstruct the new superblock Hamiltonian H'^^, which 
in turn has the same form as the original HgB in Eq-(il) substituting all the operators by 
their primed versions. This statement can be explicitly checked by considering the set of 4 
superblocks . Firstly, the new H'^^ has a contribution coming from the truncation of each 
of the 4 superblocks, each of them contributing with an A'-operator as in Eq. (|67D . Secondly, 
H'sB picks up two more contributions coming from the horizontal and vertical interaction 
between superblocks, which we denote by and H^. Thus, in the CBRG- method H'^^ is 
renormalized as, 



H'sB 



A' 



\ 



A' 



A' 



(single superblock contribution) 



(^. 



+ 



+ 



^LR J 



+ 



B' 



C'r 



DU 



B' 



D 



I 



V 



and altogether we arrive at the previously stated result of Eq.(|65l). 



(70) 
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Similarly we may proceed with the renormalized interblock Hamiltonians H ^^^J^ (|68D and 

H'^sB^Ib (HD up with the same form for them as the original ones. 

This ends the implementation of the CBRG-method for the 2D-tight-binding model. 

In Table 5 we collect our CBRG results for the first 4 lowest lying states for a chain of 
A^ = 4x4x4^ = 65536 sites. Comparison with the exact results gives a good agreement. 
We have also data from truncations with blocks of = 9 and rtg = 16 sites which enforce 
this statement. Moreover, notice that the first excited state is a doublet as in the exact 
solution. 

Another important result of our CBRG-method is that the {n1 + n\)/N'^ scaling law 
for the energy of the (ni, n2)-excited states of a square lattice of length N is reproduced 
correctly. In fact, from data of the Ug = 16 sites truncation for the first-excited-state energy 
we can extract the corresponding 1/A^^-scaling law which turns out to be, 

E(™)(iV) = cgi^G]^, c^CBRG = 9.T365, (N oo) D=2 Free BC's (71) 

while the exact value of the proportionality constant c is Cexact = tt^ = 9.86. This amounts 
to a 1.3 % error. 

Likewise, we may obtain the full {nl + n\)/N'^ scaling law for the whole set of 15 excited 
states and we find, 

Et!,nS\N) = cgLc^^^^l^, ccBHG = 7-9074, (AT oo) D=2 Free BC's (72) 

which now amounts to a 10.5 % error. 

As in the ID Free- Free case, we can determine critical scaling exponent 9 (0). For a 
truncation scheme 16 — 4 we find, 

e = 1.99999981 D=2 (73) 

which clearly supports the scaling laws introduced above. Notice again (see Table 5) that 
our CBRG method gives the exact (within machine precision) energy of the ground state. 
This is true for every step of the RG, as was proved in for ID. 
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We can also perform the wave function reconstruction of the excited states in the two- 
dimensional real space. This is achieved by a two-dimensional extension of the reconstruction 
equation p9|). As an illustration of how the CBRG method performs with this matter, see 



22 1 . The qualitative real-space form of the excited-state wave functions are captured by the 



CBRG procedure. 

With this discussion we close the first part of these notes which have been devoted to new 
develoments of the Real-Space RG method revolving around the new ideas brought about 
by the Density Matrix RG method. 



VIII. STANDARD BRG FOR THE ID AF HEISENBERG MODEL 

In the remaining sections we shall return to the standard BRG methods to deal in an 
analytical controlled fashion with models which include many-body interactions unlike the 
free models considered in the first part of these notes. 

The arquetypical model we shall be dealing with is the Heisenberg model which will be 
studied in ladder systems and 2D lattices (square, honeycomb). Our point will be that 
even with the old-fashion BRG can be useful to retrieve the correct qualitative phyisics 
when properly implemented. To this end we shall be needing some results concerning the 
Heisenberg model in one dimension which will be basic. 

Let us recall the standard BRG method for the AF Heisenberg-Ising model whose Hamil- 
tonian is given by: 

N-l 

Hn = JY. {Sp^+, + SjSjVi + ASISJ+i) (74) 

where A > is the anisotropic parameter and J > for the antiferromagnetic case. If A = 1 
one has the AF-Heisenberg model which was solved by Bethe in 1931. If A = one has 
the XX-model which can be trivially solved using a Jordan- Wigner transformation which 
maps it onto a free fermion model. For the remaing values of A the model is also solvable 
by Bethe ansatz and it is the ID relative of the 2D statistical mechanical model known as 
the 6-vertex or XXZ-model. 
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The region A > 1 is massive with a doubly degenerate ground state in the thermodynamic 
hmit N —>■ oo characterized by the non-zero value of the staggered magnetization, 

"^st = (^E^|(-m (75) 

j 

The region < A < 1 is massless and the ground state is non-degenerate with a zero 
staggered magnetization. The phase transition between the two phases has an essential 
singularity. 

We would like next to show which of these features are captured by a real-space RG- 
analysis. The rule of thumb for the RG-approach to half-integer spin model or fermion 
model is to consider blocks with an odd number of sites. This allows in principle, although 
not necessarilly, to obtain effective Hamiltonians with the same form as the original ones. 
Choosing for (^) blocks of 3 sites we obtain the block Hamiltonian: 

= l{[Si + S2 + S,]' - (5i + S^y - 3/4} + eiS^S^, + S^SI) (76) 

e := A - 1. 

If e = the block Hamiltonian Hb is invariant under the SU{2) group and according to 
the introduction to this section, we should consider the tensor product decomposition: 

-0-(g)- = -©-©- (77) 
2 2 2 2 2 2 ^ ^ 

The particular way of writing Hb given in Eq. (^) suggests to compose first Si and S3 and 
then, the resulting spin with 5*2. The result of this of this compositions is given as follows: 

||^) = |TTT) Eb = J/2 (78) 
1^, i) = i=(| Tit) + I ITT) + I TTi)) Eb = J/2 (79) 



h'2'^ y/2 



,-)i = -=(|TTT)-UTT)) Eb = Q (80) 
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1^, ^)o = ^(2| Tit) - I ITT) - I TTi)) Eb = -J (81) 

Hence for A = we could choose the spin 1/2 irrep. with basis vectors ||, |)o aiid ||, — |)o 
in order to define the intertwiner operator Tq. 

However, if A 7^ the states (|7|) -(|T]) are not eigenstates of (|76D . The full rotation group 



is broken down to the rotation around the z-axis. The states ||, |) and ||, |)i are mixed in 



the new ground state which is given by: 

T Vl + 2x2'"'2'2'' ' '"■^'2'2 



+ ^) = \ m \. \ + V2a:|^, \)) (82) 



where 



and its energy is, 



(83) 

8 + A + 3VA2 + 8 ^ ^ 



J. 



Eb = --[A + v/A^] (84) 



along with its | — |) partner. This are now the two states retained in the RG method. To 
be more explicit, we have 

I + ^) = / / [(2^ + 2)1 Tit) + (2a: - 1)1 TTi) + (2x - 1)| ITT)] (85) 
2 J6(l + 2a;2) 



I - ^) = - / / ,J (2^ + 2)1 TTi) + (2x - 1)1 iiT) + (2x - 1)1 Tii)] (86) 
2 ^^6(1 + 2x2) 

The intertwiner operator Tq reads then, 

To = | + ^)('Tl + |-^)ai (87) 

where | T)' and | |)' form a basis for the space V = C^. The RG-equations for the spin 
operators Si {i = 1,3) are then given by 

T^S^To = CS'- 2 = 1,3. (88) 
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T^S!To = eS'- z = l,3. (89) 

T^S^To = CS'- ^ = 1,3. (90) 

where etc are the renormahzation factors which depend upon the anisotropy parameter 

by, 

_ 2(l + a;)(l-2a:) 
^ ^ • 3(1 + 2x2) ^ 

^ -"3(1 + 2x2) ^^^^ 
Observe the symmetry between the sites i = 1 and 3 which is a consequence of the even 
parity of the states ( ^5]) -(^). 

The renormahzed Hamiltonian can be easily obtained using Eqs.(|88D-(p2D and (|74D , and 
apart from and additive constant it has the same form as H, namely, 

T^HnU A)To = jCbU A) + H^/sif, A') (93) 

where 

J' = in' J (94) 

A' = (p^A (95) 

Iterating these equations we generate a family of Hamiltonians if^3m(^*-^), A*^™)). The 
energy density of the ground state of Hn in the limit — > cxd is then given by, 

rp CO -| 

hm ^ = e^«^ = V -i-eB( J(-), A(-)) (96) 

where initially J^^-* = J, A*^°^ = A and Eqs.(0) -(^) provide the flow of the coupling 
constants. 

The analysis of Eq. (|9^) shows that there are 3 fixed points corresponding to the values 
A = (isotropic XX-model), A = 1 (isotropic Heisenberg model) and A = oo (Ising model). 
The properties of these fixed points are given in Table 6. 
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The computation of e^'^ in this case is facihtated by the fact that (0) becomes a geometric 
series at the fixed point. The exact resuhs concerning the models A = and A = 1 are 



extracted from references ||28| and ^y^. The case with A ^ cxo is exact because the states 
I ± |) given in ( pSf ) - ( p6D tend in that hmit to the exact ground state | tit) and | HI) of 
the Ising modeL As a matter of fact, 

1 + ^) ^A^oo I TiT)-^! TTD-^I ITT) 

I - ^) =^A^oo -| iti) - ^1 lit) - ^1 Tii) 

The region < A < 1 which flows under the RG-transformation to the XX-model is 
massless since both J(") and A^") go to zero. We showed at the begining of this section that 
all this region is critical (a line of fixed points) and therefore massless. The RG-equations 
(0) - (|95D are not able to detect this criticality except at the point A = 0. Only the 
masslessness property is detected. 

The region A > 1 which flows to the Ising model is massive and this follows from the 
fact that the product J(™)A(™) goes in the limit m — oo to a constant quantity j(°°)A*^°°^ 
which can be computed from Eqs. (|93) -(p5D and ([87|), 



t(oo) A (oo) _ TT 4 (1 + XmY 

"iio 9 (1 + 2x^)2 ^^'^ 

where Xm is given by ( ^3]) with A replaced A^'^\ This quantity gives essentially the mass 
gap above the ground state and also the end-to-end or LRO order (Long Range Order) given 
by the expectation value ■ ^(A^))! in the limit N oo. 

In summary, the properties of the Heisenberg-Ising model are qualitatively and quanti- 
tatively well described in the massive region A > 1 while in the massless region < A < 1 
one predicts the massless spectrum but no criticality at each value of A. This latter fact 
is rather subtle and elusive. One would like to construct a RG-formalism such that the 
Hamiltonian Hn{A) would be a fixed point Hamiltonian for every value of A in the range 
from -1 to 1. 
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The phase transition between the two regimes is correctly predicted to happen at the 
value A = 1. This is a consequence of the rotational symmetry, namely at A = 1 the 
system is SU (2) invariant and the RG transformation has been defined as to preserve this 
symmetry. When A 7^ 1 the SU{2) symmetry is broken and this is reflected later on in the 
RG-flow of the coupling constant A. 



IX. RG FOR HEISENBERG SPIN LADDERS 

Much of the El Escorial Summer School has been devoted to the nowdays very active 
field known as ladders systems (spin, t-J, Hubbard ...), see ||3^, @], 0. What does the 



BRG method have to say on these systems? We again emphasize that this is a technically 
simple method which produces qualitative correct results when properly applied. Later it 
is possible to look for numerical accuracy using DMRG, second order RG (see appendix) or 
some other means. 

The Hamiltonian of a Heisenberg spin ladder with ni legs, each of length is given by, 

'Hladder 'Hleg ~\~ T~trung 

ni N 

'^ie9 = J2 E ^S,(n)-S,(n+1) 

a=l n=l 
ni-1 N 

^™n5 = E E ^' S,(n) ■ Sain + 1) (98) 

a=l n=l 

where Sa{n) are spin-1/2 matrices acting on the a-th leg at the position n = 1, . . . , N, and 
J is the intraleg coupling constant while J' is the interleg exchange coupling constants, both 
beign positive to guarantee AF spin ladders. 

We shall concentrate on the uniform Heisenberg ladders with no staggering. There are 
many examples that can be worked out but to be concrete we shall pick up the 3-leg ladder 



system |^0[. In advance, what we are going to obtain is the RG-flow towards the strong 
coupling limit of spin ladders. This is an alternative to the determination of the RG-fow 
using bosonization as performed by H.J. Schulz |Q. To apply the BRG we need to set 



up what is the block Hamiltonian Hb which we do by forming blocks of 3 sites each along 
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every leg and located one block on top of another as in fig. 7 In this fashion we are selecting 
a subset of couplings from the whole spin ladder Hamiltonian in (pH]). The remaining terms 
involving links with neighbouring blocks make up for the interblock Hamiltonian Hbb- We 
shall not write down explicitely the analytical expressions for Hb and Hbb as it is quite 
clear what is meant simply by looking at fig. 7. 

Now it is aparent that the standard results of the previous section are at work for spin 
ladders. Notice that Hb is made up of small block Hamiltonians of 3 sites as in (^) whose 



eigenstates and energies are already computed in (|78|)-(pl|). The renormalization process 
goes through all the way by truncating the block states to the lowest eigenstates, i.e., the 
spin doublet |i,±i)o, with energy cq = —J- In this case the embedding operator T^") for 
each block is nothing but the projector Pq"-* onto these states; denoting 

iV/3 

Pb=T{ Pr (99) 



a=l 



then, to first order in J the renormalized or effective Hamiltonian acting on the states left 
out after the truncation is simply, 

H,^ = Pb{Hb + Hbb)Pb (100) 



Using the embedding operator (pTl) we arrive at 



PbHbPb = +yeo = --J (101) 

In this case the renormalization of the block Hamiltonian gives the identity because the two 
states retained within each block are degenerate by rotational invariance. The renormaliza- 
tion of the interblock Hamiltonian is also simple if we observe that Hbb contains products 
of spin operators belonging to different blocks, say S" ■ where a ^ j3 denote neighbouring 
blocks and i,j = 1,2,3 are the intrablock labels used in Eq. (pSf). Then, according to Eq. 
(p9|) we have, 

P^Sr ■ S?Pb = PB{Pt^StPt\Pl^^^^^Pl^^^)PB (102) 
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Hence we only need to know how the spin operators renormahze within each block onto the 
new spin operators. By symmetry arguments, the renormalization spin factor denoted by 
must be the same for the 3 components of the spin operators, i.e., 

(p(«)sapW)=^^S', (103) 

where denotes the spin 1/2 operator acting on the effective spin 1/2 subspace of the 
a*^-block. The renormalization spin factors are known from (pll) -(p2l) to be given by, 

ei = 6 = ^ 

6 = -I (104) 

Now we are ready to compute the renormalization of the interblock Hamiltonian Hbb- 
According to Eqs. (|102|) and ( p.03| ), the renormalization of the horizontal couplings between 
blocks of Hbb (see Fig. 7) is given by, 

J S^t^ ■ Sf ) J ^1^3 S,(n')S„(n' + l) = ^J Sa{n')Sa{n' + 1) (105) 

while the 3 vertical couplings between two blocks are renormalized to 

J' (Si") ■ Sf ) + S^") ■ Sf + si") • Sf )) J' + + il) S(") • 8(^5) 

= J'S,(n')-S„+i(n' + l) (106) 

We have then obtained that the renormalized Hamiltonian ( |100| ), apart from the constant 
term (|101|) , is the as the original ladder Hamiltonian, but with length A^/3 and the following 



renormalization of the coupling constants. 



J' — > J' (107) 



Hence the ratio J' / J increases as, 

J' 9 J' 



J 4 J 
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after each step of the RG showing that J' / J = oo is a stable fixed point which controls the 
behaviour for all values of J and J', while J' / J = is an unstable fixed point. Had we 
chosen blocks made up of more than 3 sites we would have obtained essentially the same 
result. The RG method for the simple spin chain (i.e. ni = 1) where first obtained in 
reference [^. According to Eq. ( |107D if we start in the weak coupling regime J' /J ^ 1, 



after sufficient iterations of the RG we would get an effective Hamiltonian in the strong 
coupling regime where we can apply the arguments of section II to derive the nature of the 
low lying spectrum of the theory. 

X. REAL-SPACE RG APPROACH TO THE QUANTUM 2D-AF HEISENBERG 

MODEL 

In this section we present a real-space RG treatment of the quantum two-dimensional 
Heisenberg antiferromagnetic model with arbitrary spin S. Most of the work using real-space 
RG methods has been devoted to one-dimensional problems. This is very useful because 
it is crucial to have an aproximate method which gives good results for both ID and 2D 
problems, for it is known that mean field theory methods fail in low dimensional problems. 
In this regard we have recently shown that the use of quantum groups in combination of 
real-space methods in ID captures the essential features exhibited by the exact solutions 



of models such as Heisenberg and ITF [T^, ||T6[. Nevertheless, the main reason which has 



prevented the applications of the real-space RG in 2D quantum lattice Hamiltonians is the 
rapid growth of the number of states to be kept in a reasonable scheme of truncation of 
states in dimensions higher than one. 

We shall be using the Block RG method in our study of the 2D Heisenberg model. This 
version of the RG method is suitable to achieve fully analytical treatments of interacting 
many-body problems. The reason for searching for complete analytical approaches as op- 
posed to purely numerical studies relies on the necessity of having a qualitative understand- 
ing of the mechanisms responsible for the different behaviors exhibited by the Heisenberg 
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model itself and for its connections to more complicated related Hamiltonians such as t-J 
and Hubbard where the understanding of the doping effects is a big issue at stake. In order 
for there to be a completely analytical RG treatment in 2D we need a juidicious choice of 
the states to be kept as we shall see pl| . 

Despite of some initial controversies there is by now sufficient theoretical and experi- 
mental evidence for the existence of antiferromagnetic long range order (AF LRO) in the 
2d spin 1/2 Heisenberg antiferromagnet (and references therein). This property has 
been observed in parent compounds of hight-Tc materials such as La2Cu04 [^. From a 



theoretical point of view this means that the strong quantum fluctuations implied by the 
low dimensionality and low spin do not destroy completely the Neel order, as it happens 
| 36| in Id. Though there is no a satisfactory physical explanation of this fact, which may be 



important regarding the interplay between antiferromagnetism and superconductivity upon 



doping. The RVB scenario originally proposed by Anderson while yielding an ap- 

pealing picture of the ground state, does not explain the presence of AF LRO. This type of 
order may however be incorporated a posteriori in long range RVB ansatzs of factorized form 



with predictions similar to the ones obtained using Quantum Monte Carlo methods |^ 



and variational plus Lanczos techniques [^. A class of physical systems where the RVB 



approach may be actually realized is in spin ladders with an even number of chains |^ , ^ 
The previous works leave still room to investigate in more depth the interplay between the 
RVB scenario, or more generally "valence bond scenarios" , and the AF order present in the 
2d AF-magnets, described by the AF Hamiltonian H = JJ2{i,j) Sj ■ Sj. 

We have proposed a new scenario where the valence bonds, instead of resonating as in the 
RVB scenario, rotate around their ends under the influence of the AF background. To test 
this idea we propose a variational ground state in which the bonds rotate but do not resonate 
among themselves. We shall start by considering how the quantum fluctuations affect the 
classical Neel state. This is also the starting point of the spin wave theory (SW), which we 
would like to use for comparison of our theory. An important ingredient of our construction 
is the use of real-space RG techniques, which allows us to obtain exact analytical results for 
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any value of the spin S of the model (S is integer or half-integer and in the discussion above 
S=l/2.) The advantage of using a real-space RG method is that one can treat in an exact 
manner the local quantum fluctuations of the classical Neel state. By this we mean that 
if we divide the square lattice into blocks of 5 sites each, as in Figs. 8,9 and 10, then the 
Heisenberg Hamiltonian restricted to the blocks can be solved exactly. The ground state for 
every block is a spin 3/2 irrep. (if S=l/2) which is obtained by forming a singlet (bond) 
between the spin at the center and the ones surrounding the center (Fig. 8) According to the 
RG method, the spin 3/2 can be chosen as an effective spin for the renormalized lattice which 
now has N/5 sites. We shall show later on that the interaction between those effective spins 
3/2 (or 3S more generally) is also governed by an AF-Heisenberg model with a renormalized 
coupling constant. Hence the RG procedure can be iterated yielding a series of effective 
spins which ultimately goes to its classical value, i.e. infinity! {S 3S ^ 3'^S oo.) 
We thus obtain in an economical and simple way the important result that the 2D AF- 
Heisenberg models belongs to the universality class of the 2D classical Heisenberg model. 
For a sigma model derivation of this result see ref. PS[. The rotating bond picture puts in 



correspondence various approaches to the 2D AF-Heisenberg model. 

Let us begin our approach by considering the cluster of 5 spins 1/2 of Fig. 8 a). The 
configuration showed in Fig. 8 a) is the exact ground state of the Ising piece of the Heisenberg 
Hamiltonian, given by //^ = JX]j=i,...,4 '^qS'/ , where Sq and are the third component of 
the spin operators at the center and the i*^ position off the center respectively. As soon as 
the "transverse" Hamiltonian H^y = JJ2i=i,...,4{SQSi+SQSf) is switched on, the down-spin 
in the middle starts to move around the cluster, and a valence bond between the center and 
the remaining sites is formed in a s-wave (/ = 0) symmetric state as shown in Fig.8 b). Other 
rotational states with / 7^ may appear corresponding to excitations (/ being the orbital 
angula momentum of the bond). An alternative description of this state is given by first 
combining the 4 spins sourrounding the center into a spin 2 irrep, which in turn is combined 
with the spin 1/2 at the center yielding a spin 3/2 irrep with energy eo = —3 J/2. If instead 
of the spin 1/2 at each site there is a spin S the previous analysis can be easily generalized 
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as follows: the ground state of the AFH Hamiltonian of the 5-cluster has total spin 3S and is 
obtained by first combining all the surrounding spins into a spin 4S, which in turn becomes 
a spin 3S after multiplication with the spin S at the center. In a certain sense this state 
can be viewed as the formation of bonds between the center and its four neighbours. After 
applying several steps of the real-space RG, as we shall see below, new bonds are generated 
between sites at longer distances apart. Thus our valence-bond scenario is a type of long 
range valence bond state. 

To study the AFH model in the entire square lattice we begin by first tesselating this 
plane using the cluster of Fig. 8 as the fundamental cell (see Fig. 9). Notice that the centers of 
the 5-cluster form a new square lattice with lattice spacing a' = ^/5a. Given this tesselation 
we can apply the standard RG method of replacing clusters of spins by an effective spin 
13| , p^ . This method has been applied for the Id AFH model by Rabin pO[ for clusters or 



blocks with 3 sites, obtaining a ground state energy with an error of 12%. The effective spin 
of every 3-block in Id has spin 1/2. In our case, as we have discussed above, the effective 
spin of the 5-blocks have spin 3S and the energy per block equal to cq = —JS{4:S + 1). The 
effective spins 5" = 35* interact by means of an effective Hamitonian which to first order 
in perturbation theory can be derived if we know the renormalization of the spin operators 
Sa ^aS',a = 0, 1, ... ,4. 

The renormalization spin factor can be shown to be given by the sum = 

^ Emo,mi,...,m4"^"(Cmo,mi,...,m4)^ subjcct to the Constraint ELo"^a = 35. C^^^rm,... is 
the CG coefficient which describes the ground state of spin 3S in terms of the 5 original 
spins S, whose expression is a product of 4 standard CG coefficients. The C,a satisfy the sum 
rule I]a=o^a = 1- We arrive at the following result, 

1 65 + l[(2g)!p 
^"^^ 3^8S + 1[(8S)!]2 

^,tm, "m(5-m,)!(5 + m,)![-2^ + EfH! ^''^ 
where if a = then mo = 3S — X^i'^i- K follows that the renormalization factors for the 
four external spins in the 5-block are all equal = ^2 = ^3 = ^4 = ^{S), while that of 
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the central spin is determined by the sum rule. Amazingly enough the sum ( |109|) can be 
performed in a close manner yielding, 

as) ^ I |i| (110) 

For spin S = ^ one obtains ^(|) = ^. Moreover, Eq. ( |110D correctly reproduces the classical 
limit lim5'_,oo ^("S") = | (recall S = S"^'^ = |S" = C,ciS'). Notice also that the value for S' = | 
is already close to the classical value. 

The RG-equations for the spin operators Si i = 1,2,3,4 allows us to compute the renor- 
malized Hamiltonian H' which turns out to be of the same form as the original AFH Hamil- 
tonian. In fact, we arrive at the following RG-equations, 

TV 

H'(N, S, J) = -JSUS + l) — 

5 

+H{f,3S,3^'^{S)J) (Ilia) 



N' = —, S' = 3S, J' = 3^^(S)J (111b) 
5 



where the first contribution in Eq. ( |llla| ) comes from the energy of the blocks. As 3C,'^{S) < 



1, the flow equation ( |1 lib] ) implies that the coupling constant flows to zero J^'^^ g -^j^j^j^ 



means that the AFH model remains massless for arbitrary value of the spin S. This fact 
allows us to compute the density of energy eoo(5') (per site) as the following series, 

oo 

eoo(5) = -| E ^J^'''^S^^\A X 3"^ + 1) (112a) 

n=0 

^(n+1) _ 3^(n)^ j(n+l) _ 3^2 ^^HyW ^^^2b) 



Using eqs. ( |112a| ) and ( |112b| ) we can compute the ground state energy of our variational 



RG state for any value of the spin S. In particular for S=l/2 we get the value e^o = -0.5464. 
This value has to be compared with the "exact" numerical result -0.6692, which is obtained 



using Green-function Monte Carlo methods PDI, and the spin wave value which is -0.6703. 



The difference between our result and the Green Function MC or SW is quite big and around 
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0.12. To clarify the origin of this departure we have considered the semiclassical expansion 
of Eq. ( |112a| ) and compare it with the standard formula of Anderson and Kubo ||49|| , 



e^^ = -2^(^ + ^ + ...) (113a) 



e^^ = -2S{S + 0.158 + ^^^ + ■■■) (113b) 



oo 



The important observation is that the term linear in S is absent in our formula ( |113a| ). The 
reason for this is that we are using a first order RG method for which the ground state 
energy follows from the formula Eqs = (\E'o|-f^|^o); where |\I^o) is the variational ground 
state constructed by the RG method. Now it is easy to see that taking |\l/o) to be simply the 
Neel state one has to go to second order perturbation theory (PT) to get a linear term in 
S, which turns out to be given by 3/4 + 1/32 + 0{1/S). It is clear that the "missing energy 
" 0.12 is due to this peculiarity of the first order PT. To remedy this one should implement 
the RG method with second order PT. In ID and for S=l/2 this can be done, obtaining for 
the ground state energy density e^o = —0.4530 which is comparable in precision with the 
spin wave result -0.4647 (recall that the exact value is -0.4431.) The latter computation in 
2D is much more involved but it is expected to yield a result close to the spin wave result. 

In order to have a better insight into the physics of the model it is convenient to compute 
the staggered magnetization M = {j^Y.j{~'^ySj)- We have been able to obtain a closed 
formula for arbitrary spin S which is capable of analytical study. To this purpose, we use 
the RG-equality for V.E.V. {ipolOltpo) = {iP'q\0'\iP'q) for renormalized observables O' in the 
ground state and divide the sum in M into 5-block contributions. With the help of the 
renormalization spin factors we arrive at the RG-equation for the staggered magnetization, 

Mm{S) = M:^^M^/5(3^) (114) 

The explicit knowledge of C,{S) ( |109| ) allows us to solve this RG-equation for the staggered 
magnetization in the thermodynamic limit ^ oo. In fact, as we know by now that 
the Hamiltonian renormalizes to its classical limit, we have lim5^oo ^("S") = 5*. Defining 
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M(S) = Sf{S), Eq. ( P^ ) amounts to solving the equation f{S) = ^^f{3S) subject to 
the boundary condition /(oo) = 1. Thus, we obtain the following formula for the staggered 
magnetization for arbitrary spin, 

ms) = sii -^1— (115) 

„=0 ^ -r 3-3 

This is a nice formula in several regards. For spin S* = | we get M(|) = 0.373 to be 
compared with the most accurate Quantum Monte Carlo reslult |^B| which is 0.3074 (earlier 



numerical results were obtained with Green function QMC methods and Variational 
Monte Carlo plus Lanczos algorithm [^). Other approximate methods employed so far 
lead to values of M(|) such as, e.g., spin wave theory plus l/S'-expansion to order S""^ gives 
^ 0.3069 (earlier SW results were provided by Anderson and Kubo [|^), spin wave 



theory plus perturbation theory gives 0.313, etc. Our value is close to the one found USTj 
with pertubation theory around the Ising model to order 4 which is 0.371. We can equally 
get values of the staggered magnetization for arbitrary spin. For spin S=l, our formula ( |115| ) 
gives 0.8454 to be compared with 0.8043 using SW to order ^ and 0.8039 obtained 
by Wheihong et al. p2[ using series expansions. 



Another interesting feature of our formula ( [115| ) is that it allows us to make a 1/5- 
expansion yielding the result, 

M^^{S) = S-0.2 + 0.06^ + 0{1/S^) (116a) 

M'^{S) = S -0.198 + 0(1/ S^) (116b) 

Observe the excellent agreement between the order S'^ term in both formulas. Recall that 
equation ( |116a| ) is derived using first order PT. We expect that a second order RG would 



further lower the value of M{S), in agreement with the numerical result. 

In summary we can claim that the rotating-valence-bond scenario gives a consistent and 
suggestive picture of the ground state of the 2D AF-Heisenberg model: the quantum fluctu- 
ations of the Neel state consist of rotating bonds which appear at all scales corresponding 
to effective spins which renormalize towards the classical value. 
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XI. APPENDIX: SECOND ORDER FORMALISM FOR THE STANDARD BRG 

METHOD 



The modern fashion to include correlations between blocks in the real-space RG is the 
DMRG [0]. Nevertheless, there is an old way to include those correlations which we believ 
has been overlooked in the past. It amounts to include a second order contribution to the 
BRG. Recall that the standard BRG in previous sections is a first order method from the 
point of view of Perturbation Theory (P.T.), i.e., the interblock Hamiltonian Hbb is treated 
perturbatively in first order. We can extend this treatment |^ to second order P.T. in the 
usual fashion and thus arrive to an effective Hamiltonian given by: 



where Eb is the ground state energy of the block and Pb denotes the projector onto te 
ground state of the block. This is nothing but the intertwiner operator of section 2. 

As a matter of illustration, we shall work out this formalism for the isotropic AF Heisen- 
berg model in ID using the 3-site block BRG explained in section 8. Moreover, we restrict 
to spin 1/2. Denote each block with an index a. Thus, the Hilbert space for each block Ti" 
is decomposed into 1/2, 1/2 and 3/2-spin subspaces. 



H,ff = Pb[Hb + Hbb + Hbb{1 - Pb) 



1 



(1 - Pb)Hbb]Pb 



(117) 



Eb — Hb 



^ 7^(")(i/2) © 7^!") (1/2) © nt\?>/2) 



(118) 



Correspondingly, we introduce 3 projectors onto those subspaces. 



(119) 



They satisfy the following properties that we will be useful. 



N' 



(120a) 



a=l 



N' N' 

^-PB = i-i[ptwm^-pt^) 



(120b) 
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p{a) p(/3) ^ p{0) p{a) 
m m' m' m 



(120c) 
(120d) 



where N' = N/3 and m indicates the site in each block. 

With the help of these properties, the second order contribution to H^ff in (|117|) , denoted 

(2) 

by -f^e/j, can be given the following form containing 3 types of terms: 

H Tf 

1 



H, 



(2) 



N' 



{[(p("-^)sr-^)p("-^)) ■ (Po^"^sS")pS)][(p£sr^p(")) ■ (p("+^)sr')p("+^)) 
-[(^o^^^s^pS) ■ (p("-^^)sr^)Po("+^))][(Po^"-^)sr^)po("-^)) ■ (pssr^p("))]}p^ 

N' ^ 

^ 2eo-e -e 

[(n^^^s^pS) ■ (pi"^'^sr'^i^ir/,^)][(i^ssr)p(")) ■ (p£ti)sr^)p("+^))]PB (121) 



In order to work out this expression ( |121|) towards a manageable result, we need to perform 
a renormalization of the spin operators both in both subspaces of spin-1/2 (recall that in 
sect. 8 we did it only for the lowest energy spin 1/2.) 
Let us introduce the following notation for the 4 states of spin 1/2: 

\m,(3) with m = ±1/2,/? = 0,1 (122) 

Denote by S' the effective spin-1/2 coming out of the block renormalization. Then, 



(m,/3|Si|m',/3') = (m|S'|m')(pi)/3,^/ 
with the p-matrices given by. 



(123) 





^ 2/3 -l/V^\ 




^-1/3 0^ 




' 2/3 1/v^^ 


Pi = 


^-1/V3 j 


, P2 = 


V 


, P3 = 


^1/^3 ^ 



(124) 



Thus, the spin renormalization that we were searching for is summarized in 
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(125) 



Namely, 



PoSiPo = P0S3P0 = = ^)S' (126a) 



PoSiPi = PiSiPo = = ^)S' (126b) 



P0S3P1 = P1S3P0 = (ef = ^)S' (126c) 



Upon substitution of these expressions in ( |121| ) we are led to the renormalization of the 
Hamiltonian: 



Hfff = E [^^'^ + Ji^^o. ■ S«+i + JaSa ■ S,+2] (127) 



a=l 



<r = E[^^°^ + '^J'^S.-S.-w] (128) 

a=l 



with the following numerical values, 

rfW = -1, ^ = I = 0.44 

' (129) 

d(2) = -0.104861, ^ = ^ = 0.130247, J2 = = 0.0411523 

Altogether, we end up with the following effective Hamiltonian in which the second order 
formalism employed shows up as a nearest-neighbour coupling J2, 

N' 

+ J2 ] (130) 

a=l 

with, 

d = -1.104861, Ji = 0.574691, J2 = 0.0411523 (131) 

We can now iterate this RG procedure as usual to obtain the RG-fiow equations for the two 
coupling constants Ji and J2: 



j{m+l) _ j{m) 
J2 — 
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(132) 



a = 0.57491 b = 0.44444 c = 0.041152 



The fixed points of tfiese RG-eqs. are simply, 



Upon iteration tlie system flows towards the smallest fixed point, 

= 0.076084 (134) 

This happens to be an understimation of the numerical value. 

Finally, we get a series expressing the ground state energy to second order in RG, 

oo -j^ 



eoo=ET^[-7i^^^+72Jr^] (135) 

m=0 '-' 



with 71 = -1.104861, 72 = 0.25. 

The above sum can be computed exactly by introducing generating functions Ji{x) = 
E'^=ox"'jt\ i = 1,2 and using the RG-eqs. (|T3^), 

Mx) = -, -7—,, ■h{x) = - ^^-^ (136) 

1 — ax -\- bcx"^ \ — ax + ocx"^ 



Thus 



1 71-72 



c 



31-f + f 



(137) 



and substituting the values of 71 and 72, we get 

e(^^^) = -0.453002 (138) 

This is to be compared with the exact value e^^*^*^*) = —0.4431 which amounts to a 2.2 % 
error. This results improves even the spin wave result e^*"-* = — S*^ — 0.36338s — 0.033011 = 
-0.4647 which is a 4 % off the exact value. Recall that e^^-^^^ = -0.391304 (11.6 % .) 
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FIGURES 

FIG. 1. a) Block decomposition of a square lattice into 4-site blocks. b)Squematic truncation 
of states in the BRG method associated to the previous lattice decomposition. 

FIG. 2. a) Ground state tpQ and first excited state ipi for the Hamiltonian Hpixed with fixed 
EC's, b) Ground state V'o and first excited state ipi for the Hamiltonian Hpree with free EC's. 

FIG. 3. Euilding blocks of the 3-site ERG for the tight-binding model in ID with free EC's, a) 
Ground state, b) First excited estate, c) Second excited state. 

FIG. 4. Pictorical decomposition of a given Hamiltonian H into uncorrelated A-operators, 
correlation B^- 5i^-operators and interaction C-operators according to the CERG method. Bi is 
a superblock made up of two Li and R\ blocks. 

FIG. 5. The /N'^-lmi for the first 5 excited states of the ID Tight-Einding Model for a chain 
of A/" = 12 X 2"* sites with Free-Free EC's. This is a lnE„-lnm plot. 

FIG. 6. The wave function reconstruction for the first 5 excited states of the ID Tight-Einding 
Model for a chain of A/" = 12 x 2^ = 768 sites with Free-Free EC's. We have scaled up the exact 
results by a factor of 1.23 for clarity. 

FIG. 7. Block decomposition associated to the standard ERG method applied to the uniform 
3-leg Heisenberg ladder. 

FIG. 8. a) The antiferromagnetic 5-block state, b) Formation of a rotating-valence-bond state 
upon applying the H^y part of the Hamiltonian to the AF 5-block state. 

FIG. 9. Artistic tesselation of the square lattice with 5-site blocks. 

FIG. 10. The two-dimensional square lattice tesselated by the 5-block. Dahed lines are near- 
est-neighbours in the renormalized lattice. 
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TABLES 



Energies 


Exact 


CBRG 


DMRG 


Eq 





1.1340 X 10-^^ 


1.0 X lO-'^ 


El 


1.6733 X 10-5 


1.9752 X 10-5 


1.6733 X 10-5 


E2 


6.6932 X 10-5 


7.6552 X 10-5 


6.6932 X 10-5 




1.5060 X 10-"^ 


1.8041 X 10-5 


1.5060 X 10-^ 


El 


2.6772 X 10^1 


2.9G81 X IQ- i 


2.6772 X IQ-'' 


E5 


4.1831 X 10-^ 


5.1078 X lO-'' 


4.1831 X 10-^ 



TABLE I. Exact and CBRG Values of Low Lying States for the ID Tight-Binding Model for 
a chain of TV = 12 x 2^ = 768 sites with Free-Free BC's. DMRG values are also given. 



m 


N=12 2"^ 


^{exact) ^j^^^ 


Ef^^^\N) 







12 


6.8148 X 10-2 


6.8148 X 10-2 


6.8148 X 10-2 


1 


24 


1.7110 X 10-2 


1.7375 X 10-2 


1.7110 X 10-2 


2 


48 


4.2826 X 10-^ 


4.4694 X 10-3 


4.2826 X 10-3 


3 


96 


1.0708 X 10-^ 


1.1515 X 10-3 


1.0708 X 10-3 


4 


192 


2.6772 X 10-^ 


2.9681 X 10-^ 


2.6772 X 10-^ 


5 


384 


6.6932 X 10-5 


7.6552 X 10-5 


6.6932 X 10-5 


6 


768 


1.6733 X 10-5 


1.9752 X 10-5 


1.6733 X 10-5 




» 1 




9.8080/A2 


9.8696/A2 



TABLE IL Exact and new CBRG values of the first excited state for the ID Tight-Binding 
Model with Free-Free BC's. DMRG values are also given. 
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Ti^ripr crips 


Exact 


Standard BRG 


CBRG 


TP 


i. ( ( 04 X iU 


i.O M i X iU 


i.o4uy X iu 


T? 

h,\ 


1.0U4O X iU 


4.ZD(y X lU 


i.DDOO X iU 


E2 


4.1761 X 10""^ 


4.2794 X 10-2 


4.6408 X 10-^ 


E3 


8.1831 X 10"^ 


4.3053 X 10-2 


9.1450 X 10-^ 




1.3520 X 10-^ 


4.3173 X 10-2 


1.5179 X 10-3 


E5 


2.0196 X 10-3 


4.4288 X 10-2 


2.2852 X 10-3 



TABLE III. Exact, Standard RG and CBRG Values of Low Lying States for the ID 



Tight-Binding Model for a chain of iV = 12 x 2^ = 384 sites with Free-Fixed BC's. 
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Ti^ripr crips 


Exact 


Standard BRG 


CBRG 


DMRG 


TP 


D.DOoO X iU 


O.oliD X iU 


(.Uo4o X iU 


D. ( X iU 


T? 

h,\ 


Z.DDoo X iU 


c 01 c c w 1 f\—2 
O.olOO X iU 


z.y4Uo X lu 


/.DO X iU 


E2 


5.9924 X 10-"^ 


5.8268 X 10-2 


6.3690 X 10-^ 


5.99 X 10-"^ 


E3 


1.0653 X 10-3 


5.8470 X 10-2 


1.2289 X 10-3 


1.065 X 10-3 




1.6644 X 10-3 


5.8717 X 10-2 


1.7707 X 10-3 


1.664 X 10-3 


E5 


2.3966 X 10-3 


5.9106 X 10-2 


2.7311 X 10-3 


2.397 X 10-3 



TABLE IV. Exact, Standard KG and CBRG Values of Low Lying States for the ID 
Tight-Binding Model for a chain of iV = 12 x 2^ = 384 sites with Fixed-Fixed BC's. DMRG 

values are also given. 



Energies 


Exact 


CBRG 


Eq 





9.6114 X 10-35 


El 


1.5056 X 10-^ 


1.9390 X 10-^ 


E2 


1.5056 X 10-^ 


1.9390 X 10-^ 


E3 


3.0012 X 10-^ 


3.8781 X 10-^ 



TABLE V. Exact and CBRG Values of Low Lying States for the 2D Tight-Binding Model for 
a lattice ofA^ = 4x4x46 = 65536 sites with Free BC's. 



TABLE VI. Fixed Points of the Anisotropic AF-Heisenberg Model 



A 





1 




00 


r>BRG 
^00 

pSxact 

gexact X iUU 


-0.2828 
-0.3183 
11% 


-0.3913 
-0.4431 
12% 
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Figure 7 



